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ANALYSIS OF TAILORED- INTERFACE OPERATION OF SHOCK TUBES 


WITH HELIUM-DRIVEN PLANETARY GASES 
By William J. Loubsky and John 0. Reller, Jr. 
Ames Research Center 


SUMMARY 


The perfect -gas relations which determine the conditions for tailored 
operation of a shock tube are summarized, and a computer method for calcula- 
ting real-gas effects is presented. The analysis incorporates the following 
assumptions: (a) the driver gas behaves as a thermally and calorically per- 
fect gas; (b) the test gas is in thermodynamic equilibrium; and (c) viscous 
effects are neglected. Equations are derived that enable the determination of 
the time, T, during which the reservoir properties remain constant (e.g., the 
time between the arrival of the incident shock wave and the arrival of the 
reflected expansion head at the end of the driven tube) . It is shown that a 
maximum value for t exists for a range of driven -to -driver tube length 
ratios; an optimum ratio is defined for minimum driver length at maximum r. 

A minimum value for T also exists if tailoring is to be achieved and cor- 
responds to an allowable maximum component length ratio. Charts are presented 
which facilitate the determination of the tailored conditions for a constant 
area shock tube having a helium drive and air, C0 2 , N 2 , or A as test gases. 
Driver conditions are varied to include pressures from 10 to 1,500 atm and tem- 
peratures from 300° to 10,000° K. Substantial deviations from the perfect-gas 
results are found for tailored Mach numbers greater than 8. For a given 
reservoir or ref lected-shock enthalpy, real-gas effects may reduce driver tem- 
perature requirements by 30 to 50 percent. Driver energy is reduced also, 
since real -gas effects enable the use of a shorter driver length. Steady -state 
times t, however, can be substantially less than the perfect -gas values. 
Examples of the use of the charts are included. 


INTRODUCTION 


One factor that limits the test time available in a shock tunnel is the 
arrival at the nozzle entrance of the reflected driver expansion wave. One 
method of making this disturbance negligibly small (as illustrated in refs. 1 
and 2) is to use a large volume driver section. This is not practical when 
energy requirements for the driver gas are so high that the driver gas must be 
arc -heated; the present state of the art makes it necessary for driver volumes 
to be quite small. 

A second limitation on test time in a shock tunnel is the reflection of 
pressure waves from the interface of the driver gas and test gas. 



Hertzberg et al. (ref. 3 ) were one of the first groups to report a method of 
eliminating these disturbances. The method was called the tailored -interface 
mode of operation and it was observed that test times per foot of driven tube 
should be approximately six times those for the reflected, nontailored condi- 
tion. The technique consists in selecting an initial loading pressure in the 
driven tube for a given driver condition so that no downstream wave is gen- 
erated when the reflected shock wave interacts with the driver gas interface. 
Ford and Glass reported (ref. 4) the basic interaction as early as 1954. 
Interface tailoring has since been studied by many groups. Flagg (ref. 5) 
provides a detailed analysis for perfect gases in which equations for the tai- 
lored shock Mach number are derived and results are presented for Mach numbers 
up to 8. 

The purpose of the present report is to show the extent to which equilib- 
rium real-gas considerations influence the tailored-interface condition for 
air and other test gases, as well as to provide information for general use in 
the preliminary design of facility components and experiments. A computer 
method is presented which extends the perfect -gas analysis to include real- 
gas effects when the test gas is in equilibrium. Charts are provided to 
enable the determination of the tailored-interface mode of operation for a 
constant area shock tube employing a helium driver and air, C0 2 , and A as 

test gases. The range of driver temperatures considered is from 300° to 
10,000° K and pressures from 10 atm to 1,500 atm; thus cold, resistance heated, 
and arc-heated drivers are included. From these charts and the gasdynamic 
charts in references 6 through 12, the shock tube flow and thermodynamic prop- 
erties can be determined. Equations are derived which enable the determina- 
tion of the available steady-state time in the reservoir as a function of the 
driver and driven tube lengths. The optimum and maximum component length 
ratios are derived. All results are presented In graphical form. The results 
indicate substantial deviations from the perfect -gas analysis. 

In practice, various mechanisms for energy transfer, such as viscosity, 
radiation, and nonideal wave interactions, will reduce performance below the 
results presented herein, but their consideration is beyond the scope of the 
present report. 


SYMBOLS 


a sound speed 

e internal energy per unit mass 

h enthalpy per unit mass 

L length 

M shock Mach number 

m molecular weight 
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p pressure 

R gas constant of cold gas 
T temperature 

t time 

U flow speed relative to shock wave 

u flow speed relative to test facility 
V shock wave speed relative to test facility 
x coordinate length from diaphragm station 

7 ratio of specific heats 

p density 

r steady -state time in the test gas reservoir 

Subscripts 

a condition upstream of standing normal shock wave 

b condition downstream of standing normal shock wave 

D driver section 

DT driven tube section 

L limiting value 

max maximum value 

min minimum value 

opt optimum value 

R reflected shock wave 

RHe reflected shock wave in helium 

RR reflected shock wave from interface 

S incident shock wave 



ST tailoring value of incident shock wave 

TS reservoir test gas 

-l condition ahead of incident shock wave in quiescent gas 

2 condition behind incident shock wave 

3 condition behind incident interface 

4 condition in driver section prior to diaphragm breakage 

5 condition behind reflected shock wave in the test gas 

6 condition behind reflected shock wave in the driver gas when the 

interface is tailored 

7 condition behind reflected shock wave in the driver gas when the 

interface is not tailored 

0 condition behind the reflected shock wave from the interface when the 
interface is not tailored 


Reference Constants 


p Q 1 atm = 1.013x10 s N/m 2 

R 0 286.8 m 2 /sec 2 °K for air; 296.9 n^/sec 2 °K for N 2 ; 189.2 m 2 / sec 2 °K for 
C 0 2 ; 207.8 m 2 /sec 2 °K for argon; 2078 m 2 /sec 2 °K for He 

T 0 273 - 16 ° K 

Ti 300° K 

p 1.290 kg/m 3 for air; 1.246 kg/m 3 for N 2 ; 1 - 9^0 kg/m 3 for C 0 2 ; I.780 kg/m 3 
for argon 


ANALYSIS 


The tailored-interface condition exists, by definition, when no unsteady 
wave disturbance propagates into the shock tunnel reservoir as a result of 
the interaction of the reflected shock wave and the driver gas interface; 
that is, the interface is brought to rest by the reflected shock wave. The 
general flow process is illustrated in figure 1 with a physical schematic and 
distance-time diagram for a constant area shock tube. The interface is 
"tailored" when the pressure ratio P 8 /p 5 - 1-0. 
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Perfect Gas Tailoring 


In this report a thermally and calorically perfect gas is referred to 
simply as a perfect gas. The interaction of a shock wave with a contact sur- 
face (interface) is reported for a perfect gas in reference k, and the crite- 
rion for tailoring is expressed in terms of the specific internal energy ratio 


7 4 + 1 + £5 
f3 = 74-1 Pg 
e 2 7 \ + 1 + P5 
7i - 1 P 2 


( 1 ) 


Flagg in reference 5 gives the tailoring requirement in terms of the tempera- 
ture ratio T 3 /T 2 by the expression 


£3 

T 2 


(7 4 — l)(l + 


m 1 ( 7 1 - l)(l + 


\ 


74+I P 5 ^l 

7 4 -i p 2 y 

n + 1 £ 5 ^ 

7 1 - 1 P 2 / 


( 2 ) 


From equation (2), the tailoring requirement is satisfied by a unique tempera- 
ture ratio across the incident interface for a specified gas combination and 
reflected shock pressure ratio. 


The tailored shock Mach number ^ as illustrated in reference 5, for a 
shock tube of constant area can be found in terms of the diaphragm speed of 
sound ratio by a simultaneous solution of equation ( 2 ) and the following 
relations governing the flow process 


T, 


- r 


L\ 7 1 + 1 


(Mg 2 - 1) + Ms £ 


7i + 




V 


(Mg 2 - 1) + 1 


J Ms c 





\ ( 3 ) 


p ( 37 ! - 1)M S 2 - 2( 7l - 1) 

5 _ __ _ . 

p 2 ( 7l - 1)M S 2 + 2 


For purposes of this report the strong shock limit of the above relations 
is used. The error introduced by the approximation^ which enables the tai- 
lored shock Mach number to be expressed in explicit terms ^ is less than a few 
percent for Mgip > 6 . For the strong shock approximation the tailored shock 
Mach number can be expressed as 
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( 7 X + l) 2 (7 4 - 1) - [ 274 ( 7 ! - l)(r 4 + l)(7 x + l) 3 ] l/2 
(7 4 ~ !) 2 (7! + 1) - 2y 4 (y 1 - l)(7 4 + 1) 

4 

and the tailoring requirement, equation (2), becomes simply 

T 3 = (7 4 + l)m4 
T 2 (y x + l)mi 



w 


( 5 ) 


From equation ( 4 ) it is apparent that the tailored Mach number is only a func- 
tion of the gas combination and the diaphragm speed of sound ratio . 1 When 
real -gas effects are considered, however, since the specific heats and degree 
of dissociation are functions of two thermodynamic properties, the tailored 
conditions in general as well as M ST are also dependent upon the loading 
pressure Pi- 


Real-Gas Tailoring 

General considerations .- Real -gas effects such as variable specific heats, 
dissociation, and ionization preclude explicit relations, such as equations (l) 
to ( 5 ) • la general, one must solve for the thermodynamic and flow properties 
downstream of the incident and reflected shock waves, and after the interaction 
of the reflected shock wave with the driver gas interface. In the present 
analysis thermodynamic equilibrium is assumed, and the thermodynamic gas 
models of Bailey and Lomax of the Ames Research Center are used (e.g., see 
ref. 13 for carbon dioxide). Each species is assumed to behave as an ideal 
gas. The species considered in each model are the following: (l) air: 

0 2 , N, 0, NO, N + , 0 + , N0 + , N ++ , e”; (2) nitrogen: N, N + , N ++ , e"; ( 3 ) 

carbon dioxide: C0 2 , CO , 0 2 , C, 0, C+, C ++ , 0 + , 0 ++ , e~; (4) argon: A, A + , 
A ++ , e". 

The flow properties across the incident and reflected shock waves are 
governed by equations of conservation of mass, momentum, and energy and the 
equation of state: 

Pa 4 »a U a 2 = P b + I 

»a + | U a 2 = \ + | V [ 

Pb = °b (p b’ h b ; J 


x Note that equation ( 4 ) differs from the general expression for shock 
tube Mach number by virtue of the added restraints imposed by the tailored 
condition . 


6 



The explicit form of the equation of state is determined from the assumed 
gas models . 


The driver gas is assumed to he perfect and expanded isentropically by an 
unsteady expansion wave, the governing equation being 


^3 = 


2a 4 

?4 - 1 


1 





( 7 ) 


The validity of the perfect-gas assumption for a helium driver gas is based 
upon the fact that both electronic excitation and ionization are sufficiently 
small for the driver conditions investigated, to have a negligible effect on 
the thermodynamic properties. The results of reference l4 illustrate this 
point . 

Computational procedure .- An IBM 7094 Fortran II computer program is used 
in the analysis, the input information being the conditions in the driver and 
driven tube sections before the diaphragm breaks . Equations (6) and ( 7 ) are 
solved simultaneously for the conditions behind the incident shock wave and in 
the fully expanded driver gas; the pressure and velocity are matched across 
the incident interface. The reflected shock properties are then found in a 
similar manner, with the condition that u 5 = 0 . The input loading pressure 
p l is selected sufficiently low that the reflected wave from the driver-gas 
interface is a shock wave. Equations (6) are solved simultaneously for the 
conditions behind the transmitted shock wave in the driver gas and the 
reflected shock wave in the test gas, using the interface boundary conditions 
u 7 = u Q and p y - P Q • The pressure ratio P q /p 5 i s then formed and a new 

loading pressure p is determined from the equation p = (pVpJp, 

1 mew Q ^ 1 initial 

The entire procedure, beginning with the incident shock wave, is repeated 
until the pressure ratio P s /P 5 < 1-005- When this criterion occurs, the 
interface is considered to be tailored. The various forms in which the con- 
servation equations are expressed in the procedure as well as the methods of 
iteration are presented in appendix A. 


Component Lengths and Steady-State Time Relations 

The purpose of tailored- inter face operation is to achieve a maximum 
steady-state time, t, in the test- gas reservoir. It can be shown that t is 
a function of driver and driven tube lengths and that an optimum component 
length ratio (L dt /L;d) t exists for any given initial conditions in the 

driver and driven tube. This component length ratio is an optimum in the 
sense that it is the only ratio for which both the driver energy required is a 
minimum and t is a maximum. This is best illustrated in the following 
sketches: 
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Sketch (a) 


From sketch (a) it is seen that 
(L OT /L£)) 0 pt occurs when the reflected 
shock wave, the reflected expansion 
headland the expansion tail all inter- 
sect at the same point on a distance- 
time diagram. When the component 
length ratio is less than the optimum 
value (sketch (b)), t is limited by 
the disturbance generated from the 
interaction of the reflected shock 
wave and the driver gas expansion tail. 
Hence for a fixed driven tube length . 



increasing the driver length so that 
(L dt /L d ) < (L DT /L D ) opt requires more 

driver energy (larger volume) with no 
increase in the steady-state time, t. 
For a fixed driver length (and energy) , 
(L dt /L d ) becomes less than (Lp^Lp)^^ 

when the driven tube length is too 
short. The result is that t < tu 


max 


and the test-slug length is also 
reduced. Hence the situation depicted 
in sketch (a) is the optimum condition. 


For component length ratios equal 
to or greater than the optimum value 
Sketch (b) the reservoir steady-state time is 

limited by the arrival of the reflected 
expansion head from the driver. Increasing the component length ratio beyond 
the optimum value decreases t until the reflected expansion head intersects 
the interface before the reflected shock wave, thus interfering with the tail- 
oring process. Hence, for tailoring to occur the component length ratio must 
he less than a maximum value (L DT /Lp) max (see sketch (c)). The allowable 
(bpT/^B^max therefore implies a minimum value, T m: j_ n . The dependence of t 
on the component length ratio is shown qualitatively in sketch (d) for fixed 



Sketch (c) 



Sketch (d) 
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thermodynamic conditions in the driver and driven tube. It is seen that the 
curve for t depends on which component length is allowed to vary. From the 
consideration of slopes on a distance-time diagram, relations for the variables 
above can be derived in terms of flow, wave, and sound speeds in the various 
regions, and the driver-gas specific heat ratio. These relations are derived 
in appendix B and are written in the following form: 


(Note: 7 = 7 4 in eqs . (8) through (15)) 


7+1 

2(7-1) 


L DT ] 

■ L D /opt 



a4 


Ldt 


2(u 3 + V R )(l - 


7 - 1 % 


a. 


2(7-1) 


'D /max 


1 + 


Vr 

V ST. 


a. 


(8) 


(9) 


r 



(10) 


(in 
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+ 1 + 


'RHe 

a 6 


2 a 4 - (7 - 1 ) 


u 


a 4 + ” p - 7 a 3 + V RHe 


1 - 


7-1 U 3 
2 a 4 


7+1 

>(7-1) 


3,4. 


'( 12 ) 


valid for 


ADT\ < ^DT < Adt/ 

V^D /opt V J D /max 


•where 


V 

RHe 




(13) 


and 



( 1*0 


1 /p 

a 4 = (7R4T4) 


( 15 ) 
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Equation (12) enables the determination of the available time t for any 
component length ratio between the optimum and maximum values. Equations ( 8 ) 
through ( 12 ) are valid when the test gas is treated as a perfect gas as well 
as when real-gas effects are considered. The variables as, U 3 , P 5 /p 2 > an( ^- V R 
are found in the charts of references 6 through 12 . 


The length of the test-gas slug contained between the stationary inter- 
face and the driven tube end-wall for the tailored condition is determined 
from the relation 


or in terms of velocities 


■ J TS 


= d^DT 


J TS 


(^ST ~~ 

(% + u 2 )V ST 


"DT 


(16) 


RESULTS AND DISCUSSION 


Comparison of Perfect- and Real -Gas Calculations 

The variation of tailored shock Mach number with driver temperature and 
initial loading pressure in the driven tube is shown in figure 2 for the sev- 
eral test gases considered. The effects of initial pressure are seen to be 
small. Comparisons with the perfect-gas strong-shock approximations illus- 
trate a substantial deviation in tailored Mach numbers greater than Mg T = 8 ; 
for example, for T 4 = 8000° K, p x = 0.001 atm in air, the real-gas Mach number 
is 23 .I as compared to 19 *^ for the perfect-gas value, an increase of 19 per- 
cent. In figure 2(d), real-gas effects in argon result in a lower Mgrp than 
the perfect-gas value for 1000° K < T 4 < k000° K. This behavior is due to the 
rapid onset of ionization in the test gas behind the reflected shock wave with 
increasing temperature. The relatively sudden transition to the ionized state 
reflects the existence of only sparsely populated excited electronic states 
that is characteristic of this gas (see ref. 10) . The results for reflected 
shock pressures shown in reference 10 are consistent with present findings. 

It should be noted that in figure 2 the driver pressure p 4 is not arbitrary 
and must be determined from figure 3* However, since the effects of p^^ are 

small, then the effects of p on the tailored shock Mach number are also 

4 

small . 

Figure 3 illustrates the effect of driver pressure and temperature on the 
initial loading pressure that is required in the driven tube to achieve inter- 
face tailoring. It is seen that p is more sensitive to the driver condi- 
tions at lower driver temperatures. For C0 2 (fig* 3( c ))j the curves for 
driver temperatures from 3000 ° to 6000 ° K lie close together while those for 
4000° and 5000° K are seen to cross. Although not completely understood, this 
behavior is attributed to relatively large changes in chemical composition of 
the test gas with temperature in this interval and to the pressure dependence 
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of this change. For this range of driver temperatures the dissociation of C0 2 
molecules is completed behind the incident shock wave and the further break- 
down of the gas into C and 0 atoms is accomplished in substantial and vary- 
ing amounts behind the reflected shock wave, as shown by the results of 
reference 9 - 

The test-gas reservoir or reflected shock enthalpy and pressure for 
tailored operation are shown as a function of driver temperature in figure 4. 
Again the driver pressure p 4 is not arbitrary and must be determined from 
figure 3 . The enthalpy is in the dimensionless form hs/RT 0 and the pressure 
is ratioed to the initial pressure p x - Again the real-gas effects are 
obviously substantial. For air, with T 4 = 8000° K, p x = 0.001 atm, the res- 
ervoir enthalpy is 37 percent higher than the perfect-gas value and the res- 
ervoir pressure is 270 percent higher. This, of course, is a direct result 
of the real -gas tailoring Mach number being greater than the perfect-gas 
value. When the reservoir enthalpy is a specified quantity, say for design 
purposes hs/RT 0 = 700 in air, then if the real-gas effects are considered 
the driver temperature required will be less than the required perfect-gas 
value, in this case by 25 percent. Similar results are observed for N 2 and 
C0 2 - With argon, however, the advantageous effects of real -gas properties are 
less pronounced. In fact, for hs/RT 0 < 500 real- gas effects are disadvanta- 
geous, reflecting again the presence of substantial amounts of ionization in 
the test gas at relatively low values of T 4 . 

The optimum and allowable maximum component length ratios are shown in 
figure 5 as a function of driver temperature with p x as a parameter. The 
real-gas length ratios are greater by as much as a factor of 2 and a surpris- 
ingly large influence of test -gas chemistry can be seen for C0 2 and argon, at 
the T 4 values noted earlier in figures 2 and 3 . From comparison with the 
perfect-gas values it is noted that, for a given dr iven tube length , consider- 
ing real -gas effects enables the use of a shorter driver; therefore, less 
total driver energy is required for a specified driver pressure and tempera- 
ture. Conversely, for a given driver length the optimum driven tube length 
must be substantially greater than predicted by a perfect-gas analysis to 
achieve a maximum real -gas, steady-state reservoir time. The intermediate 
situations where Lpij/Lp is less than optimum should be avoided since they 
have been shown earlier to have adverse effects on steady-state times and 
driver energy requirements (see sketches (b) and (d)). 

Figure 6 illustrates the effects of driver temperature and initial pres- 
sure p x on the steady-state time parameter a-jr/lp. The parameter, and 
therefore the available steady-state time t , decreases with increasing tem- 
perature and depends slightly on the initial pressure in the driven tube . 

For a given driver temperature and initial pressure, the consideration of real- 
gas effects yields a value for the parameter less than the perfect-gas value. 
The deviations are less pronounced for the optimum values of LdtAt) than 
for the maximum values. For example, at T 4 = 8000° K, p x = 0.001 atm from 
figure 6 (a) the deviation in the parameter is only 8 percent for the optimum 
case but is 40 percent for the maximum. This is so for two reasons: First, 

in the optimum case the controlling feature is the velocity of the primary 
expansion tail. While this disturbance moves much faster relative to the main 
shock wave in the real-gas flow than in the perfect gas, the added driven tube 
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length required for real-gas optimization is compensatory and t is not 
greatly different. Second, in the case of maximum Lp<p/Lp the length of the 
test-gas slug (see eq. (l6)) is dominant and is sufficiently reduced in the 
real gas that even the greater driven tube length does not overcome the dif- 
ference, and t is substantially less- In summary, if driver length is 
fixed, the increase in driven tube length to meet the real-gas optimum condi- 
tion largely compensates for the time reduction (per foot of tube) due to 
real-gas effects. If, however, the driven tube is fixed and the driver is 
shortened to minimize energy requirements for real -gas operation, then the 
combined results of figures 5(a) and 6(a) reveal the steady-state time 
T opt ^ 0.40 T opt • F° r a long driver, when Lpip/Lp is less than 
real perfect 

(L DT /LD) opt time does not increase so long as the driven tube length is 
fixed, as was indicated in the Analysis section. This last condition might 
not be restrictive if the relatively slow change of pressure in the primary 
expansion tail can be tolerated (see sketch (b)). 


Use of Charts 

The use of the charts will be illustrated by the following examples . The 
first example is to determine the conditions required for tailored operation 
of a shock tunnel for a given set of driver conditions, say T 4 = 4000° K and 
p 4 = 400 atm, in a helium-air gas combination. From figures 2(a) and 3(a) the 
required shock Mach number is Mg>p = 14-7^ and the initial driven tube loading 
pressure p 1 = 0.10 atm. The reservoir or reflected shock enthalpy and pres- 
sure are found to be I 15 /RT 0 = 365 or h 5 = 28-72X10 6 joules/kg, and 
p 5 /pi = 3350 or p 5 = 335 atm, from figure 4(a) . These values as well as the 
remaining reservoir properties and the conditions downstream of the incident 
shock wave can be found from the charts of references 6, 7, and 12. Pertinent 
values for this example are u 2 = 4.62 mm/psec, Vp = 0.843 mm/psec, and 
as = 2.24 mm/psec. From equations (8) through (ll) or from figures 5(a) and 
6(a), (L DT /L D ) opt = 2.5, ( L DT/ L D) max = 7*3> ( a i T /L D ) mx = 0.152, and 
(a 1 T/Lj 3 ) n] q ri = 0.062. Therefore for a driver length, Lp, of 3-0 m, 

(■^DT^ opt — 4*5 m, (LpT) max ~ 21.9 ^"max ~ l*3l4 msec, and T ]nin — 8.536 msec. 
The length of the reservoir test slug is found from equation (l6); 

( L Ts) opt = °- 111 m > ( L Ts)max = °' 3 28 m- 

The second example is to determine the conditions required for tailored 
operation of a shock tunnel for a He-C0 2 gas combination for a given reservoir 
enthalpy and pressure, say hs = 28.3OXIO 6 joules/kg and p 5 = 9° atm. From 
figure 4(c) the driver temperature is T 4 = 3000° K and the loading pressure 
is p 1 = 0.01 atm, found by a trial and error selection of p x in the figure. 
From figures 2(a) and 3(a) the required shock Mach number is Mgrp = 19-0 and 
the required driver pressure is p 4 = l4o atm. The same procedure as in the 
first example is used to determine the remaining reservoir properties, con- 
ditions behind the incident shock wave, component lengths and steady-state 
times t. 
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CONCLUDING REMARKS 


An analytical and numerical study has been made of shock tube performance 
in the tailored -interface mode of operation, with primary emphasis on the 
influence of equilibrium real -gas effects in the test medium. The effects of 
viscosity, energy radiation, and nonideal wave interactions have not been con- 
sidered in this study, and to this extent the specific results will be at 
variance with experimental findings. 

The perfect-gas relations for the tailoring requirement are summarized 
and a computer method is presented that extends the perfect-gas results to 
include equilibrium real-gas effects in the test gas. Equations are derived 
for determining the time during which the reservoir properties remain con- 
stant. Charts are presented that facilitate the determination of the flow and 
thermodynamic properties for tailored operation for a constant -area shock tube 
employing a helium driver and air, C0 2 , N 2 j or argon as test gases. Driver 
conditions are varied from 300° K < T 4 < 10,000° K and 10 atm < P 4 < 1,500 atm. 

The results deviate substantially from the perfect-gas results when real- 
gas effects are considered, especially for tailored Mach numbers greater than 
8 . With respect to driver length and energy requirements these deviations are 
very favorable. For a given driver pressure and temperature the real-gas 
tailored Mach number, and hence the reservoir enthalpy, is considerably higher 
than the perfect-gas value. This result is important in an experiment (e.g., 
heat transfer) for which the reservoir enthalpy level is a given quantity. Con- 
sidering real-gas effects can reduce the driver temperature requirements by as 
much as 30-50 percent, with a corresponding reduction in driver energy require- 
ments. In addition to the reduction in driver temperature requirements for a 
given reservoir enthalpy, a further reduction in total driver energy is appar- 
ent from the fact that a shorter driver is required for a given driven-tube 
length . 

It is shown that a maximum steady-state time, t , exists and corresponds 
to the optimum component length ratio. A minimum steady-state time also 
exists, if tailoring is to be achieved, and corresponds to the maximum allow- 
able component length ratio. The time, t, during which the reservoir proper- 
ties remain constant is reduced from the ideal -gas value. The consideration 
of real-gas effects is a necessity when conditions are determined for tailored- 
interface operation at Mach numbers above 8 and driver gas temperatures in 
excess of 1000° K. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., March l 8 , 1966 
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APPENDIX A 


COMPUTATIONAL PROCEDURE AND ITERATION METHODS 


An IBM 709^ Fortran II program incorporates as input information the 
following thermodynamic properties in the driver and driven tube: p Ti, 7 

Ri^ p 4 ^ T 4 , 7 4 ^ and R 4 . The driver gas speed of sound and limiting velocity 
are calculated from 

a 4 = (7 4 R 4 T 4 ) 1/2 (Al) 


and 



(A2) 


Equation (A2) is necessary in the iteration method for matching pressure and 
velocity across the incident interface. 


'////////////// ////////////,. 

v s u r° 

'////////// //HT// /'7T7T/T77, 

Sketch (e) 


Pa 


h-2 



INITIAL MOVING SHOCK WAVE 

The flow geometry is illustrated 
in sketch (e) . A Galilean transforma- 
tion is used to rewrite equations ( 6 ) 
as follows: 


Pi 


+ 




-\ 



(A3) 


P a = P 2 (P 2 ' h2) 




The choice of an initial value of ua/Vg , for example, ug/Vg = 1, leads 
to first approximations for p~ an( 3- ba • From thermodynamic tables a p 2 and 
thus an improved value for u 2 /Vg is found. The cycle is repeated -until 
successive values of u^/Vg differ by a negligible amount. 
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ITERATION FOR EQUAL PRESSURES ACROSS INTERFACE 


The driver gas is expanded unsteadily from stagnation conditions to a 
pressure p . From equation ( 7 ) 


2 7 4 



(A4) 


where is given and a 4 is known from equation (A2) . The iteration con- 

sists of initially choosing two values for u 2 = U 3 and comparing the result- 
ing p 3 from equation (A4) with p 2 from the initial moving shock 
relations (A3) • The method is illustrated in sketch (f) . The curve for p 3 



Sketch (f) 


is arbitrarily fitted with a parabola through the two chosen values for u^ 
and the limiting velocity, u^* The p 2 curve is approximated arbitrarily 
with a straight line. The intersection of the two curves defines a new esti- 
mate for the velocity u^ which is used to compute new p 2 and p 3 . The 
iteration proceeds in this manner until the following criterion for conver- 
gence is satisfied: 



< 0.00005 


u 2 

777777777777777777777777777 / / 

Sketch (g) 



REELECTED SHOCK WAVE 


The flow geometry is illustrated 
in sketch (g) . Equations ( 6 ) of the 
text, after a Galilean transformation, 
are rewritten in the form 


1 6 



p 


5 


*2 




h 


5 


h-2 






p 5 = p 5 (p 5 ^ 5 ) 


(A5) 


where the properties in region 2 are known. The choice of an initial value of 
u^/Vr (e.g., = l) leads to first approximations for p 5 and h 5 . From 

the thermodynamic tables a p 5 and thus an improved value for u^/Vr is 
found. The process is repeated until successive values of u^/Vr differ by a 
negligible amount. 


REFLECTION AND REFRACTION OF THE REFLECTED SHOCK WAVE 

AT THE INTERFACE 


The flow geometry is illustrated in sketch (h) . For the refracted shock 
wave, which transmits into the driver gas (considered to be a perfect gas), 
the shock equations become 


and 



Sketch (h) 


Before passing 
through interface 


After passing 
through interface 


P 7 - P 3 + P 3 ( u 3 + VRHe) 2 (l " p^) 

^7 = u 3 + V RHe 
P 3 


(A6a) 

(A6b) 
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(a6c) 


1a_± 1 £t _ i 
£z = ^4 ~ I P 3 
P 3 ?4 I _ .£7 
7 4 - 1 P 3 


The reflected wave from the interface can be a shock wave, a Mach wave, or an 
expansion wave, depending upon the initial loading pressure in the driven tube. 
When it is a shock wave, the governing equations are 


p 8 P 5 + 1 


h 8 = h 5 + iffi 


1 - 
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£5 = x £e 
P 8 ' " Vl 


RR 


(at) 


P 

a 


°e (p e’ h s ) 




The wave reflected from the interface will be a shock when u y = u 8 > 0. 
The boundary conditions are p 7 = p and u = u Q • The procedure for solving 
equations (A 6 ) and (A7) is as follows: 


(1) Let u 8 = 0, as a first guess. 

(2) Assume a P 7 /p 3 and calculate a P 7 /p 3 from equation (A 6 c) . 

( 3 ) Calculate Vpg e from eqiation (A 6 a) . 

(4) Calculate P 7 /p 3 from equation (A 6 b) and repeat steps (2) to (4) 
until the successive values of P y /P 3 differ by a negligible amount. 

(5) If the input p^^ is sufficiently small, the resulting p is 
greater than p • The computer now solves equations (A 7 ) for the guessed 
value u Q and compares p with p . If p 4 p a new estimate for u 
xs made . 


( 6 ) Using the new u 8 , repeat steps (2) to ( 5 ) until successive values of 

u 8 differ by a negligible amount, thus determining the conditions in regions 

(7) and ( 8 ) . The result is that p =p,u = u_ > 0, and p /p >1. 

7 8 7 8 8 5 


DETERMINING THE INITIAL DRIVEN TUBE LOADING PRESSURE FOR 

TAILORED OPERATION 


For a given set of driver conditions, the in put. p 2 is selected suffi- 
ciently small that a shock wave is reflected from the interface as illustrated 
in sketch (h) . New values of p x are successively computed from 
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p 


x new 



CAS) 


and the entire process beginning with the incident shock wave is repeated 
until the following criterion is established: 

1.000 < < 1.005 

As p /p lj u 7 = u s -*■ 0 and the interface is considered to he tailored. 
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APPENDIX B 


DERIVATION OF COMPONENT LENGTH AND RESERVOIR 
STEADY-STATE TIME RELATIONS 

GENERAL EQUATIONS 


The equations are derived in terms of wave and flow velocities and the 
driver specific heat ratio. Variables used in the derivations are defined in 
the distance-time diagram below. 



Sketch (i) 


The time for which the reflected expansion head from the driver wall and 
the expansion tail intersect, t 3 , is found by the method of reference 15* For 
a perfect gas it can be shown that 

(Note: 7 = 7 4 i n "the following equations) 


(at 2 ) 3 = / a^Y 

(at 2 ) 4 " W 
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which, after rearrangement, becomes 


7+1 

2(7-1) 


, 3-4 

t 3 = t4 la^ 


(Bl) 


The driver length is 


Lj) - a4t4 
or 

l d 

t4 = a? (B2) 

The equation for the isentropic unsteady expansion of a perfect gas from 
stagnation conditions can he written as 


54 = ( ± _ 7 - 1 , 
as \ 2 a 4 y 


(B3) 


Combining equations (Bl) , (B2) , and (B3) results in the intersection time 




7+i 

2(7-1) 


(b4) 


From the sketch the following equations can be written: 
Lj)t ~ I3 t ljj|_ l<i t* 1,3 


(B5) 


L3 - (u 3 - a3)t3 


L i " V RHe^ t3 * - tp _ ^(^i* - ^3*) 

= (U 3 + a3)(t3* - t 3 ) 


1*2 = V R (t ± - t 2 ) = V R t ± 




%_ L 

V ST DT 


(B6) 

(B7) 

(B8) 


(B9) 


where L DT /Vg T = t 2 
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L3 + 

From equations (B5), (B9) , and (BIO) 

Vr- 


1 + 


t-i = L 


'DT 


V ST 


(BIO) 


U3 + Vr 

\ ^ 

Equating (B 6), (B7), and (B8) to equation (BIO) yields 

(U3 - a 3 )t3 + V^jj e t3* - + (U3 + a 3 )t 3 * - (U3 + a 3 )t ; 


(Bll) 


3 ~ 


1*3^ 


From equation (B3) , the above equation can be rearranged as follows 


t 3 * = 


u 3 + V 


RHe 


^4 + ^ + Vgg 



2a 4 - (7 - l)u3 

a 4 + “ p 7 ^3 + Vjyjg 


Combining equations (B 4), (Bll) , and (B12) gives 


t 3 


(B12) 


ts* = 


U3 + V- 


RHe 


> a 4 + ~ V 7 U3 + 



2a 4 - (7 - l)ug 


3 - 7 

a 4 + — - — U3 + VRH e 


7-1 U; 

2 K 


7+1 

2(7-1) 


^4 


l d (B13) 


The time for which the reflected expansion head and the tailored- 
interface interact is found by solving equation (B7) for t-j_*, 


t, = ( 1 + 

1 V a 6 


t 3 * - 


cl <2 


t . 


(Bl 4 ) 


The time, t, during which the reservoir thermodynamic properties remain 
constant is the interval between the arrival of the reflected expansion head 
and the incident shock wave in the reservoir. From equation (B9) , 
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T = ~ t • + t^* 
a 5 1 1 



(B15) 


When equations (Bll ) , (Bl4 -) } and (B15) are combined, t becomes 


T = 


'Vr 

a 5 



l dt 


+ i + 


fes)t 3 * 

a 6 / 


(Bl6) 


The general equation for the time t in terms of known driver and driven 
tube lengths and the various wave, flow, and sound speeds is found by combin- 
ing equations (B13) and (Bl6), 


t - l d < 


V R V RHe 



+ 1 + 



U-3 + V- 


RHe 


3-7 

a 4 + U«n + V] 


2 




l dt 

l d 


+ 1 + 


^RHe 

a 6 


2a 4 - (7 - l)u 3 
3-7 

a 4 + — g U3 + Vftjjg 


- 7+1 H 

_ 7-1 2 ( 7 - 1 ) 

2 a 4 


a 4 


(B 17 ) 

Equation (B17) holds for any L ET /Lj) such that Lj_* > 0 and Lj_ > 0. To evalu- 
ate equation (B17 ) > a 6 and V^ e must be calculated from the Rankine-Hugoniot 
relations. From the transfer of velocity coordinates in sketch (j), the ideal 
Rankine-Hugoniot relation is used to calculate the reflected shock velocity in 
the driver gas , VR£ e • 



T 7777777J777777T77777t7T/// 

Moving shock 


V//////////////////////////////////A 


u 3 + V, 


RHe 


v RHe 


! V RHe 


V7777TTT7TTT7T777777TT77TTT7777T7777, 

Standing shock 

Sketch (j) 


U 2 

u 5 


^3 + v RHe 
v RHe 


1 + 


7 + 1 Pt 
7-1 Pi 


7 + 1 + £s 


7-1 
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Solving for V^jj e , 


Vi 


RHe 



-x 


(B18) 


The speed-of-sound ratio across the reflected shock in the driver gas is 

P=- N 1/2 



+ 1 , ^5 

- 1 Pp 


i + Es 


7 - 1 P 2 , 

Combining equations (B3) and (B19) yields 


a 6 - 



1/2 


(B19) 


(B20) 


l/2 

where a 4 (/RT 4 ) • The pressure ratio 

£5 _ £5 £1 
P 2 P^, P 2 

where p 5 /p p^ /p , and the remaining quantities Vggi, Vg, 1 x 3 , a 5 can be 

found from the charts in the references. 


SPECIAL CASE EQUATIONS 


Optimal Relations 


Maximum time t can be achieved when the reflected shock wave, reflected 
expansion head, and expansion tail intersect at the same point. Prom 
sketch (i) T ^y and thus L et /L d | .. occur when t 3 * = t 3 , Lj_* = 0. From 
equations (B?) and (B13) 


t 3 * = 



7-1 ^3) 
2 a 4 y 


7+1 

2(7-1) 




2 a 4 - (7 - l)ug 

a 4 + ^ 2 7 113 + V RHe 



7-1 ua\ 
2 a 4 / 


7+1 n 


8,4 


l d 


2A 


* 


♦ 


I 



which when rearranged yields the optimum component length ratio 


l dt 

Ld 


7+1 


U 3 - a 4 + 


VRHe^ * Vr)(i - 2 -j J: 2}) 


7+1 

2(7-1) 


opt 


(u<3 + VRHe)^! + &4 


(B21) 


From equations (b 4) and (Bl 6 ) with t 3 * = t 3 , the maximum time t becomes 


T max h> 


'Vr _ VRHeM 1 + V S T 

, as a s J Vus + V R 



• Lj DT 

l D 


opt 


+ 1 + 


v RHe 
a q 


— 

7+1 -i 


(1 - y ~ 1 Ha) 

H 

1 

K 

c\l 


V 2 3 - 4 , J 



a4 



- 

- 

j 


(B22) 


To evaluate equation (B22) one must assume a driver length. 

Minimal Tailoring Relations 

An upper limit on the component length ratio, is imposed since 

the reflected shock wave must interact with the interface prior to the 
reflected expansion head if tailoring is to occur. From sketch (i), the limit 
results when t 3 * = tj_* = tj_, Lj_ = 0. From equations (Bll) and (B13) the 
maximum allowable component length ratio becomes 

7+1 


l dt 

h) 


2(ua + V R )(l 


i(7 _ i) 


max 


Vr 

1 +^]a 4 

VST 


(B23) 


From equations (Bll) and (Bl 6 ) with t 3 * = t^, the minimum time t becomes 


t . = Ln 

m in & 


1 + 



1 + 


V R 

V ST 


U 3 + V* 


R 



DT 




(B2k) 


max 


Again one must assume a driver length. 
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Figure 1.- Schematic and distance-time diagrams for a constant area shock tube 
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Figure 2.- Continued. 
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Figure 2.- Continued 
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Figure 2.- Concluded. 









Helium driver pressure, p 

















Dimensionless enthalpy, (h 5 /R 0 T 0 ) x 10 



(a) Air. 

Figure k.- Reflected shock enthalpy and pressure for tailored conditions 

(P 4 not arbitrary) . 


Ratio of reflected shock to initial pressure, (p 5 /p, ) x 10 




















Component length ratio, 



(c) Carbon dioxide. 
Figure 5-- Continued. 
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Steady-state time parameter, a ( r/L 



k6 


(a) Air. 

Figure 6.- Steady-state time at tailored conditions; optimum and maximum 

Lpij/Lp ratios. 
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